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ABSTRACT 

In  this  paper  we  consider  the  problem  of  calculating  force  targets  for  a 
collection  of  n  fingers  which  grasp  a  two-dimensional  object  at  known  posi- 
tions, at  which  the  normals  to  the  surface  are  also  assumed  to  be  known  at 
least  approximately.  If  the  points  at  which  the  fingers  touch  the  body  do  not 
allow  a  positive  grip  to  be  exerted  (i.e.  a  grip  in  which  the  fingers  hold  the 
body  in  equilibrium  by  exerting  friction-free  forces  in  the  directions  of  the 
corresponding  inward-directed  normals),  it  is  appropriate  to  find  the  smallest 
coefficient  of  friction  for  which  it  is  possible  to  assign  a  set  of  finger-force 
targets  which  hold  the  object  at  equilibrium  and  such  that  each  individual 
force  lies  within  the  corresponding  cone  of  friction.  We  present  an  algorithm 
for  this  problem  which  runs  in  time  0(n  log^/i  log  log  n).  We  also  present 
another  algorithm  for  preprocessing  the  given  data  so  as  to  allow  fast  compu- 
tation of  the  desired  coefficient  of  friction  for  the  case  in  which  one  needs  to 
balance  any  given  'query'  external  force  and  torque.  Finally  we  discuss 
simpler  variants  of  our  techniques  which  are  likely  to  be  more  efficient  when 
the  problem  is  solved  for  a  small  number  of  fingers. 

1.   Introdnction 

The  geometric  problem  analyzed  in  Sections  2  and  3  of  this  paper  arises  in  the  following 
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context.  Wc  consider  a  robot  'hand'  of  some  design,  consisting  of  several  independently  mov- 
able, force  sensing  'fingers'.  For  simplicity  we  shall  suppose  that  these  fingers  move  as 
points  in  the  plane  and  that  each  can  sense  the  vector  of  force  that  it  applies  to  any  object 
which  it  contacts.  (An  experimental  planar  manipulator,  called  the  four-finger  manipulator, 
having  these  characteristics  has  been  constructed  in  our  laboratory;  see  [DLSS].)  The  fingers 
are  assumed  to  be  quite  stiff,  but  also  assumed  to  have  some  small  degree  of  compliance 
(relative  to  the  positions  to  which  they  have  been  commanded  to  move);  non-zero  compliance 
is  necessary  if  the  fingers  are  to  grip  perfectly  rigid  bodies  without  infinite  forces  arising. 
Under  these  assumptions,  what  patterns  of  control  should  we  apply  to  the  motions  of  such 
gripping  fingers,  in  order  to  cause  them  to  move  a  gripped  object  along  a  specified  path, 
while  steadily  applying  a  specified  level  of  grip  force? 

Although  this  central  pragmatic  question  defines  the  context  addressed  by  this  paper,  we 
will  not  seek  to  discuss  all  of  its  aspects  here.  Instead,  we  shall  merely  outline  the  control 
approach  developed  in  the  companion  papers  [Sc],  [DLSS],  which  is  as  follows:  It  is  assumed 
that  the  system  to  be  controlled  ('hand',  plus  the  object  which  it  is  gripping)  is  moving  slowly 
enough  for  all  inertial  effects   to  be  negligible.  Then 

A)  The  best  available  model  of  the  system  to  be  controlled  is  analyzed  using  the  principles 
of  elementary  mechanics,  leading  to  the  definition  of  a  'good'  set  of  gripping  points  and  of  a 
force  target  for  each  of  the  gripping  fingers.  This  is  the  force  which  the  finger  is  supposed  to 
exert  during  the  required  motion  (assuming  that  the  external  forces  and  torques  on  the  object 
will  remain  constant  during  that  motion). 

B)  Step  (A)  prepares  for  a  control  loop,  which  aims  to  keep  the  finger  forces  actually  meas- 
ured during  an  experiment  close  to  their  target  values.  If  the  discrepancy  between  model  and 
actual  system  is  not  too  large,  this  should  cause  the  object  O  to  follow  a  path  close  to  that 
desired,  while  the  forces  exerted  on  various  specified  points  of  O  remain  near  the  target 
values  specified  in  subphase  (A). 
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The  (finger-)  force  targets  used  in  subphase  (B)  are  derived  by  relating  them  either  to 
various  user  specified  external  force  targets,  or  to  measurements  of  external  forces  (we 
assume  that  those  measurements  can  be  performed  using  the  fingers  themselves).  The  rela- 
tionship between  externa]  and  target  force  used  during  subphase  (A)  in  any  particular  case 
must  reflect  the  nature  of  the  particular  manipulation  task  to  be  accomplished.  For  example, 
we  may  wish  to 

a)  Use  several  fingers  to  grip  a  freely  movable  object  O,  and  then  move  O  along  a  specified 
path  while  applying  grip  forces  of  user-specified  intensity. 

b)  Grip  an  object  O  constrained  by  various  hinges,  pegs,  and  slots,  and  move  it  along  the 
unique  path  available  to  it,  while  keeping  the  forces  that  oppose  the  constraints  on  O  as  small 
as  we  can. 

c)  Apply  grip  to  an  object  O  lying  on  a  table-top,  and  exert  force  sufficient  to  overcome 
gravity  (here  note  that  the  weight  and  center  of  gravity  of  O  will  generally  not  be  known  in 
advance). 

The  fingers  of  the  'hand'  are  assumed  to  grip  the  contacted  body  frictionally.  The  force 
and  torque  acting  on  the  body  at  any  moment  can  then  be  calculated  from  the  forces  acting 
between  the  body  and  all  of  the  gripping  fingers;  it  is  assumed  that  sensors  for  these  latter 
forces  are  present  in  the  fingers. 

If  we  model  the  gripping  fingers  as  points  which  can  move  in  2-space  and  make  (fric- 
tional)  contact  with  the  gripped  object,  the  only  relevant  aspects  of  the  geometric  situation 
are  the  positions  at  which  the  fingers  are  to  contact  the  object,  and  the  directions  of  the 
tangents  to  the  object  at  these  grasp  points.  Suppose  first  that  the  gripping  points  have 
already  been  chosen.  Besides  the  evident  equilibrium  conditions  (force/torque  balance),  to 
prevent  slip  we  must  insist  that  the  force  F,  exerted  by  the  «-th  gripping  finger  must  lie  in  the 
cone  of  friction  at  the  corresponding  grip  point  x,.  Let  \x.  be  the  coefficient  of  friction  (which 
we  assume  to  be  the  same  for  all  fingers),  and  let  v,*(ji)  and  v^(ji)  be  the  two  boundaries  of 
the  friction  cone  at  the  i-th  point,  so  that  v* {yL)  =  ni  +  [it,  and  vr(fi.)  =  n,-ji»,,  where  n,  (resp. 


t,)  is  the  unit  normal  (resp.  unit  tangent)  to  the  gripped  object  at  the  i-th  grip  point.  Then 
slip  will  be  avoided  as  long  as  F,  belongs  to  the  cone  spanned  by  v,*(n,)  and  v~  (ii)  for  all  i. 
The  two  equilibrium  conditions  that  constrain  F,  are 

^F,  =  0        and        '^F.^x,  =  0  . 
These  conditions  can  be  written  as 

where  ptii^)  (resp.  pTCm-))  is  the  3-dimensional  vector  (vf  (jjL),(vf  (ji,)xa:,)3  )  (resp. 
(vf  (ji),(v~(ji)xX|)3  )).  We  require  a  solution  of  the  system  (*)  satisfying  a*sO,  a,~aO  for 
all  i  (and  of  course  not  all  the  a,'s  can  be  zero);  hence  the  a,'s  can  be  normalized  to  have  a 
sum  equal  to  1,  making  it  clear  that  what  is  needed  is  to  express  the  zero  vector  as  a  convex 
combination  of  the  vectors  p^  (\l),  p~ (ii). 

Suppose  that  n  fingers  are  involved.  Then  (inequalities  aside)  the  2n  coefficients  o,*, 
a~  arc  constrained  only  by  the  three  linear  conditions  (*),  which  are  far  from  determining 
them  uniquely.  To  resolve  this  ambiguity  in  an  advantageous  manner,  we  propose  to  use  the 
following  heuristic.  The  nominal  coefficient  ja  of  friction  must  be  regarded  as  somewhat  unc- 
ertain, so  it  is  appropriate  to  look  for  the  smallest  ji*  for  which  a  positive  solution  a* ,  a.~  of 
the  system  (*)  exists,  and  then  to  use 

F,  =  arvr(jt*)  +  arv,-(tt*) 
as  the  target  force  for  the  I'-th  finger.  (Note  that  the  solution  which  best  protects  against  fric- 
tion coefficients  smaller  than  assumed  is  also  the  solution  which  best  protects  against  angular 
uncertainties  in  the  tangent  directions  at  the  grip  points.)  When  |i.  =  0,  the  friction  cone 
degenerates  to  the  vector  n,  normal  to  the  object  surface  at  the  point  x,,  and  our  problem 
becomes  that  of  expressing  the  zero  vector  as  a  convex  combination  of  the  vector 
^.=  (n„(/iiXx,)3  ). 

In  the  case  of  three  gripping  fingers  there  will  generally  exist  no  positive  a,  satisfying 
2ai^i=0  since  the  three  vectors  N,  will  generally  not  span  a  triangle  passing  through  the  ori- 
gin of  3-space.    (Such  a  solution  can  only  exist  if  the  three  lines  normal  to  the  object  surface 
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at  the  three  grip  points  are  concurrent.)  However,  if  four  fingers  are  used,  X^a,A',  =  0  may 
have  a  solution,  and  in  this  case  F,  =  a,Ni  gives  the  appropriate  force  target  for  the  »-th  finger 
(note  however  that  we  can  multiply  all  these  forces  by  a  common  factor  to  achieve  whatever 
total  grip  force  is  desired). 

In  the  preceding  discussion  we  have  assumed  that  the  gripping  positions  x,  are  fixed  in 
advance.  However,  a  more  challenging  variant  of  the  problem  is  to  find  a  good  set  of  grip- 
ping positions,  given  some  approximate  representation  of  the  object  to  be  gripped,  which  can 
be  acquired  from  visual  or  other  sensory  data.  Suppose  for  simplicity  that  the  object  boun- 
dary is  represented  as  a  simple  polygon  of  n  sides.  Let  e,  be  the  i-th  side  of  that  boundary, 
and  consider  the  vectors  v*(fjL)  =  n,-i- p./,  and  v~  (^L)  =  n,-  \Lt,,  where  n,  (resp.  t,)  is  a  unit  vec- 
tor normal  to  e,  and  pointing  into  O  (resp.  tangent  to  e,).  The  collection  of  all  forces  propor- 
tional to  one  of  these  vectors  which  are  applied  at  arbitrary  points  along  e,  generate  a  cone  of 
force/torque  vectors  in  3-space  which  is  easily  seen  to  be  generated  by  just  the  four  vectors 
Pij(i>^)>  Pi7 (►*■)>  J  -  1.2,  where  p,)(m-)  (resp.  p,7(fi.))  is  the  3-dimensional  vector 
(v/'(»».),(v/'(M.)xe,y)3  )  (resp.  (vr(M.),(vr(pi)x«,;)3  )),  where  «,i.  «.2  are  the  two  endpoints 
of  «,.  In  other  words,  we  seek  the  smallest  ii*  for  which  the  origin  can  be  expressed  as  a 
positive  (or,  actually,  convex)  combination  of  the  4n  points  ^^(ji*),  pfiit-'),  i  =  1,  .  .  .  ,n, 
j  =  1,2  in  3-space.  Once  this  p,*  has  been  found,  application  of  Caratheodory's  theorem 
implies  that  the  origin  can  also  be  expressed  as  a  convex  combination  of  just  four  of  these 
vectors.  These  vectors  will  then  define  'optimal'  gripping  positions  for  the  fingers.  See 
[MSS],  [MY]  for  a  similar  analysis  of  the  synthesis  of  good  grips. 

These  observations  define  the  problem  studied  in  the  following  sections.  In  the  sequel 
we  analyze  the  problem  assuming  a  large  number  n  of  grasping  fingers,  and  obtain  asymptot- 
ically efficient  algorithms  for  its  solution.  This  assumption  may  not  be  realistic  in  practice, 
but  it  facilitates  a  rather  non-trivial  analysis  of  the  geometric  structure  of  the  problem  and 
enables  us  to  establish  interesting  connections  between  our  problem  and  several  sophisticated 
techniques  in  computational  geometry.  Moreover,  the  analysis  in  Section  3  leads  to  an  algo- 
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rithm  which  is  quite  suitable  for  practical  use,  perhaps  with  some  minor  modifications  that 
trade  fancy  data  structures  for  cruder  ones  that  are  easier  to  use.  The  concluding  section  4 
addresses  the  issue  of  more  pragmatic  implementations  of  our  techniques. 

Before  closing  the  introduction,  we  note  that  many  variants  of  the  grasping  problem 
have  been  studied  recently.  In  addition  to  the  papers  [MSS],  [MY]  mentioned  above,  the 
problem  has  been  also  considered  in  [BFG],  [HA],  [LK],  [Ngl],  [Ng2].  However,  the  treat- 
ment of  frictional  grips  as  done  in  this  paper  is  new  to  the  best  of  our  knowledge. 

2.    An  efficient  grip-finding  algorithm  for  a  free  object 

The  first  abstract  problem  to  which  the  preceding  discussion  leads  is  as  follows:  Let 
Piiii.),  ■  .  .  ,p„{ii)  be  n  points  in  3-space,  depending  linearly  on  a  real  non-negative  parame- 
ter fi.  Let  Ar(^.)  =  conv  {p  i(|j.),  .  .  .  ,Pn(M-)}  be  their  convex  hull,  and  let  C  (ji)  be  the  cone 
generated  by  ^(ji)  with  the  origin  as  the  apex  (i.e.  C(ji)  =  {Xv  :  v  €  ^(|t),  X>0}).  We 
assume  that  as  ^l  increases  so  does  C((j.).  Our  problem  is  to  find  the  smallest  n*  for  which 
the  origin  lies  in  Ar(fi.*).  In  this  section  we  present  an  algorithm  for  this  problem  which  runs 
in  time  0(n  log^n  log  log  n). 

We  begin  with  the  following  simplification  of  the  problem:  For  each  x,y  and  each  jjl 
define 

Af'(x,y,M.)  =  minp,(ji,)-(x,y,  1) 

M'(x,y,y.)  =  minp,(|jL)(x,y, -1) 

Note  that  there  exists  a  point  (x,y)  for  which  either  M  "^  (x.y.jt)  >  0  or  Af  ~(x,y,ji.)  >  0  if  and 

only  if  the  origin  lies  outside  K(y.).  Thus  our  task  is  to  calculate  the  two  quantities  it,*  and 

ji"  which  are  respectively  defined  as  the  smallest  p.  for  which 


resp. 


max  M  ■*■  (x.y,  ji)  s  0 

•t.> 


max  M    (x,y,jji)  s  0  . 


The  larger  of  these  two  ji's  is  the  quantity  p,*  that  we  seek. 


We  only  need  consider  the  problem  of  finding  jjl*,  since  \i~  can  be  found  in  the  same 
way.  For  simplicity,  write  M*  =  M.  First  fix  »a,  and  consider  the  subproblem  of  calculating 
max  M{x,y,iL).  Our  assumptions  imply  that  if  this  maximum  is  positive  (resp.    nonpositive) 

we  have  \i.*>\i  (resp.    p,"^S(i). 

For  each  fixed  \i.  we  face  a  linear  programming  problem  in  R',  solvable  in  linear  time 
using  an  elegant  technique  of  Megiddo  (see  [Mel],  [Me2];  cf.  also  [Dy]).  However,  we  will 
need  to  combine  this  technique  with  another  technique  (also  due  to  Megiddo;  see  below)  for 
finding  the  appropriate  p.. 

The  algorithm  makes  use  of  two  subprocedures  FINDMAX(jjl)  and  LINETEST(L,^.). 
FINDMAX  either  establishes  that  the  maximum  of  the  function  M  =  M(x,y,ii.)  for  a  specific 
(L  is  positive,  or  locates  this  maximum;  LINETEST  is  called  by  FINDMAX  to  do  the  same 
for  the  restriction  of  A/  to  a  specified  line.  Both  FINDMAX  and  LINETEST  exploit  the  fact 
that,  since  M  (x,y,  p.)  is  the  minimum  of  a  finite  collection  of  functions  linear  in  x  and  in  y,  it 
is  concave  in  these  variables.  Hence  it  has  a  unique  local  maximum,  and  M  assumes  its  max- 
imum value  either  at  a  unique  point  which  is  the  intersection  of  three  or  more  planar  faces  of 
the  polyhedral  graph  of  M,  or  at  all  the  points  of  a  horizontal  edge  of  this  graph,  or  at  all  the 
points  of  a  horizontal  face  of  this  graph.  FINDMAX  will  call  LINETEST  repeatedly,  to  find 
the  maximum  of  M  over  various  specified  lines  L.  Suppose  that  this  restricted  maximum  value 
V  is  attained  at  a  point  p  ^  L.  Once  LINETEST  has  located  p  it  is  easy  to  determine  the  side 
of  L  in  which  the  global  maximum  of  M  lies.  Indeed,  if  Af  s  v  everywhere  near  p  along  the 
line  L'  orthogonal  to  L,  then  since  M  is  concave  p  is  clearly  a  local  maximum  of  M,  from 
which  it  follows  that  v  is  the  global  maximum  of  M.  Otherwise  M  exceeds  v  at  some  point  q 
near  p  along  L' .  Then  plainly  there  can  exist  no  point  q'  on  the  side  of  L  opposite  to  q  which 
is  such  that  M(q',ix,)  a  v,  since  if  there  were  such  a  q'  the  concavity  of  M  would  imply  the 
existence  of  a  point  along  L  at  which  M  had  a  value  larger  than  v.  These  facts  will  be  used 
implicitly  in  what  follows. 
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The  strategy  of  FINDMAX  is  to  cut  down  the  number  of  planes  in  3-space  which  must 
be  considered  in  calculating  the  maximum  of  M;  each  of  the  main  steps  of  the  procedure  cuts 
away  a  fixed  proportion  (approximately  1/16)  of  the  planes  remaining  at  any  given  moment. 
This  is  done  by  combining  these  planes  into  disjoint  pairs  p,  finding  the  line  L(p)  of  intersec- 
tion of  each  such  pair,  and  then  finding  pairs  p  of  planes  for  which  we  can  make  an  a  priori 
determination  of  the  side  of  L(p)  that  contains  the  poini  at  which  the  global  maximum  of  M 
is  attained.  One  plane  of  each  such  pair  can  then  be  discarded  from  further  consideration. 
The  routine  FINDMAX  uses  LINETEST(L,m.)  in  the  manner  already  explained  to  determine 
the  side  of  L  in  which  the  maximum  is  attained,  if  this  maximum  is  nonpositive. 

It  is  easiest  to  regard  the  routines  FINDMAX  and  LINETEST  as  receiving  specific 
values  of  ji  as  their  parameter.  However,  these  routines  are  called  by  a  main  routine  FIND- 
p,,  which  calculates  n^  by  simulating  a  'generic'  execution  cf  FINDMAX  without  knowing 
the  desired  value  of  n'  a  priori.  This  is  done  using  an  interesting  and  significant  idea  of 
Mcgiddo  [Me3],  which  is  explained  in  more  detail  below. 

Having  given  this  "bird's  view"  of  the  algorithm,  we  proceed  to  describe  each  of  its  rou- 
tines in  more  detail.  (Again  we  remind  the  reader  that  the  first  two  routines  FINDMAX  and 
LINETEST  are  adapted  directly  from  Megiddo's  linear  programming  algorithm.  Because  of 
the  sophisticated  manner  in  which  these  algorithms  arc  used  in  the  main  procedure  FIND-ji, 
we  spell  them  out  in  detail.) 

LINETEST:  Finding  the  maximam  of  Af  on  a  line  L 

First  consider  the  restricted  subproblem  of  calculating  the  maximum  of  M  on  a  given 
line  L  in  the  plane.  Because  of  its  relationship  to  the  other  routines  comprising  our  algorithm 
(FINDMAX,  FIND-jji),  we  require  the  procedure  LINETEST  which  finds  this  maximum  to 
have  one  of  the  following  three  outcomes  in  every  case:  (i)  It  finds  a  positive  value  for  the 
maximum  of  M  on  Z,.  In  this  case  the  procedure  FINDMAX  that  invokes  LINETEST  (again 
we  note  that  the  aim  of  FINDMAX  is  either  to  calculate  the  maximum  of  M  over  the  whole 
plane  or  to  establish  that  this  maximum  is  positive)  can  also  stop  and  report  back  to  its  caller 
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FIND-ji  that  the  required  ^."  must  be  greater  than  the  value  ji  supplied  to  FINDMAX  as 
parameter,  (ii)  LINETEST  determines  that  the  global  maximum  of  M  is  attained  on  L,  and 
that  this  maximum  is  nonpositive.  In  this  case  execution  of  FINDMAX  can  also  terminate, 
reporting  back  that  ji^  ^  ji.  (iii)  LINETEST  finds  the  (relative)  maximum  of  A/  on  L  to  be 
nonpositive,  but  not  to  coincide  with  the  global  maximum  of  M.  In  this  case,  it  will  also 
determine  on  which  side  of  L  the  global  maximum  of  M  lies. 

The  LINETEST  procedure  is  as  follows.  Without  loss  of  generality  we  can  assume  that 
L  is  the  X-axis,  so  that  our  problem  is  that  of  computing  max  M(x,0,n).  Since  the  graph  of 

X 

F{x)  =  Af  (x,0,n.)  is  the  lower  envelope  of  n  linear  functions  L],  .  .  .  ,L„,  it  must  be  a  con- 
cave function  of  x.  To  find  its  maximum  we  proceed  as  follows.  Group  the  lines  L|,  .  .  .  ,L„ 
into  disjoint  pairs  of  lines.  Find  the  point  of  intersection  of  the  two  lines  in  each  such  pair, 

and  let  these  points  be  X] x„    where  n'  =    ln/2j .    Find  the  median  xq  of  xj,  .  .  .  ,!„■ 

(in  linear  time),  and  (also  in  linear  time)  calculate  F^xq).  If  F(xq)  happens  to  be  positive,  we 
find  ourselves  in  the  above  case  (i),  i.e.  since  a  positive  value  of  M  has  been  found, 
LINETEST(L,n)  (and  the  procedure  FINDMAX(ji,)  which  will  have  called  it)  can  return  an 
indication  that  the  desired  jjl*  is  greater  than  the  ^l  supplied  to  LINETEST  as  parameter. 
Next  suppose  that  /"(xq)  s  0.  Let  x*  denote  the  point  at  which  F  attains  its  maximum. 
Plainly  if  the  line  L  at  which  F  (xq)  is  attained  has  positive  (resp.  negative)  slope,  then 
x*  >  Xq  (resp.  X*  <  xo);  if  L  is  horizontal,  we  can  take  x*  =  xq,  and  exit  from  LINETEST 
since  the  relative  maximum  of  M  on  Z,  have  been  found  (in  this  case  we  will  find  ourselves 
either  in  case  (ii)  or  in  case  (iii),  as  indicated  below). 

Suppose  for  the  moment  that  x*  i=  xq  and  then  without  loss  of  generality  that  x*  >  xq. 
It  follows  that  X,  <  X*  for  each  of  the  roughly  n/4  points  x,  lying  to  the  left  of  xq.  Thus,  of 
the  two  lines  in  the  pair  defining  x,,  it  is  the  line  with  larger  slope  that  will  be  higher  at  x*. 
Since  this  line  can  be  clearly  discarded  without  changing  the  maximal  value  F(x*)  of  the 
lower  envelope  function  F,  we  can  remove  /i/4  lines  from  our  original  set  without  changing 
the  desired  maximum.   After  this  we  can  re-run  the  same  procedure,   applying  it  to  the 


-  10  - 

remaining  3n/4  lines,  and  can  repeat  this  process  until  our  set  of  lines  becomes  small  enough 
for  a  concluding  direct  calculation  to  be  inexpensive.  This  gives  a  linear  time  overall  pro- 
cedure for  calculating  the  maximum  value  v  of  F.  If  this  maximum  is  positive  we  find  our- 
selves in  case  (i)  again  and  can  return  an  indication  that  ii.*>\i.  However,  if  v  is  nonpositive, 
we  can  examine  the  normals  to  the  collection  of  planes  along  which  F(x*)  is  attained  (usually 
we  expect  to  have  only  one  or  two  such  planes,  but  in  degenerate  situations  F(x*)  can  be 
attained  by  any  number  of  planes).  If  the  xy-projections  of  all  these  normals  lie  in  an  open 
half-plane  H  bounded  by  a  line  /,  then  M  is  easily  seen  to  increase  when  moving  from  x*  in 
the  direction  normal  to  /  and  pointing  away  from  H;  otherwise  x*  is  a  global  maximum  of  M. 
This  observation  allows  us  to  determine  whether  the  global  maximum  of  M  is  attained  at  a 
point  lying  at  the  upper  half  plane,  the  lower  half-plane,  or  on  the  x-axis.  Thus  LINETEST 
can  always  return  one  of  the  three  possible  results  described  above. 

This  completes  our  description  of  LINETEST. 

FINDMAX:  Finding  the  global  maximum  of  A/ 

The  second-level  routine  FINDMAX  must  either  find  the  global  maximum  of  M  (for 
our  given  p.),  or  show  that  this  maximum  is  positive,  in  which  case  ji^  >  ft.  This  routine  can 
be  constructed  as  follows.  Group  the  planes  /,(x,y)  =  p,(ji,)(x,y,  1)  into  disjoint  pairs,  and 
find  the  line  of  intersection  L,j  of  each  pair/,,  fj.  Megiddo's  technique  is  to  find  some  fixed 
portion  of  these  lines  for  which  LINETEST  can  find  what  side  of  the  line  Ly  contains  the 
point  at  which  the  global  maximum  of  M  is  attained.  If  this  can  be  done  we  can  discard  one 
of  the  planes  in  the  pair  (f,,fj)  because  M  will  lie  below  this  plane  at  its  global  maximum. 
This  same  process  can  then  be  repeated  for  the  remaining  collection  of  planes,  and  we  can 
iterate  until  either  some  positive  value  is  found  for  M  (as  in  case  (i)  of  LINETEST),  or  until 
the  desired  (nonpositive)  global  maximum  is  found  on  one  of  the  lines  L,^. 

We  claim  that  each  of  the  above  line-pruning  steps  can  be  performed  in  linear  time, 
yielding  a  linear  time  overall  procedure  for  calculating  the  global  maximum  of  M.  To  see 
this,  rotate  the  coordinate  system  if  necessary  to  ensure  that  at  least  half  of  the  lines  Ljj  have 
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nonnegative  slope,  while  the  other  half  of  the  lines  have  nonpositive  slope  (this  can  plainly  be 
done  in  linear  time  by  calculating  the  median  of  the  angles  formed  between  the  lines  L^  and 
the  positive  j;-axis).  Combine  the  lines  in  disjoint  pairs  L,j,  Lki,  putting  one  line  with  nonnega- 
tive slope  and  one  line  with  nonpositive  slope  into  each  pair.  Each  pair  of  lines  has  a  unique 
point  of  intersection,  unless  both  the  lines  of  the  pair  are  horizontal,  in  which  case  we  take 
them  to  determine  the  point  on  the  y  axis  which  is  half-way  between  them.  Next  calculate  (in 
linear  time)  the  median  y„  of  the  y  coordinates  of  these  points,  and  apply  the  line-testing  pro- 
cedure LINETEST  to  the  line  y  =  y„.  Unless  LINETEST  reaches  case  (i)  and  stops  (in 
which  case  we  know  that  M  has  a  positive  maximum),  it  will  indicate  that  the  maximum  of  M 
lies  in  one  of  the  half-planes  y  >  y„  or  y  <  y„.  (If  the  maximum  lies  on  y  =  y„  it  will  have 
been  found  by  LINETEST  (which  will  then  have  reached  terminal  case  (ii),  so  that  FIND- 
MAX  can  also  terminate.)  Suppose  then  that  the  maximum  does  not  lie  on  y  =  y„,  and  sup- 
pose for  definiteness  that  the  maximum  lies  in  y  >  y„.  Take  the  collection  of  all  intersection 
points  lying  on  the  other  half-plane  y  <  >«,  and  calculate  the  median  x„  of  their  x  coordi- 
nates. Apply  LINETEST  to  the  line  x  =  x„;  this  will  yield  another  half  plane,  call  it  x  >  j:„, 
in  which  the  global  maximum  of  M  is  attained.  Altogether  this  tells  us  the  quadrant  in  which 
the  maximum  is  attained;  in  the  specific  case  considered,  this  is  the  quadrant  at  >  x„,  y  >  y„. 
Consider  that  quarter  of  the  set  of  all  pairs  L,_,,  L*;  of  lines  whose  intersection  point  lies  to  the 
left  of  x„  and  below  y„.  For  these  lines  we  can  be  sure  that  the  global  maximum  lies  on  the 
upper  (rather  than  the  lower)  side  of  one  of  the  lines  Z.^,  L*/.  Hence  for  this  line  (say  L,j)  we 
know  which  of  the  functions  /,,  fj  is  going  to  be  larger  at  the  global  maximum  of  M,  and  thus 
we  can  discard  either  /,  or  fj  from  the  set  of  planes  under  consideration.  This  allows  us  to 
eliminate  1/16  of  the  planes  /;  as  explained  above,  a  linear  time  procedure  for  finding  the 
maximum  of  M  results. 

This  completes  our  description  of  FINDMAX  and  with  it  our  review  of  Megiddo's 
linear  programming  technique. 
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FIND-»i  -  Finding  the  optimal  ^t. 

Next  consider  the  problem  of  finding  the  quantity  p,^  defined  above.  Here  we  can  use 
another  elegant  technique  of  Megiddo,  which  is  to  run  the  above  procedure  FINDMAX  in  a 
'generic'  manner,  without  specifying  the  value  of  \i. 

Since  this  idea  is  general  and  of  broad  applicability,  we  digress  to  describe  it  in  suitably 
abstract  terms.  Megiddo's  technique,  which  might  be  termed  'symbolic  binary  search', 
applies  in  the  following  situations.  A  certain  function  F(ji)  of  a  real  variable  ji  is  defined  by 
an  algorithm  a(/4(ji,)),  whose  input  arguments  i4(jt)  are  one  or  more  arrays  whose  com- 
ponents depend  on  \i,  but  in  some  relatively  simple  manner,  e.g.  algebraically  but  as  polyno- 
mials or  rational  functions  of  low  degree.  We  seek  a  value  ja*  of  |i  defined  by  some  condi- 
tion on  F(jx),  e.g.  F(jji*)  =  0.  It  is  assumed  that  the  algorithm  a  can  be  adapted  so  that, 
whenever  it  returns  a  value  F(^l)  =  a{A{^l.))  not  meeting  the  condition  that  we  want  to  satisfy, 
it  can  tell  us  on  what  side  of  the  parameter  value  ^  the  target  n*  lies. 

An  obvious  idea  in  this  case  is  to  attempt  a  binary  search  for  ^.'.  For  example,  if 
values  ji;  and  p.^  for  which  ji;<p.*<ji,^  were  known  a  priori,  we  would  calculate 
F((fi;-i- ^t/,)/2)  to  determine  whether  ji*  lies  in  the  upper  or  lower  half  of  the  interval 
[)i,i,(jla],  and  then  iterate  this  bisection  procedure.  If  the  function  F(m.)  is  favorable,  e.g.  if  it 
is  a  monotone  increasing  or  decreasing  step  function,  this  procedure  will  terminate,  and  if 
there  are  k  critical  points  at  which  the  value  of  F((i.)  changes  and  if  these  critical  points  are 
distributed  reasonably  uniformly  within  [u/,)!^]  then  only  0(logk)  bisection  steps  will  be 
required,  yielding  an  algorithm  of  overall  cost  0(K  log  k),  where  K  is  the  cost  of  a  single 
invocation  of  the  underlying  algorithm  a.  The  problem,  however,  is  that  in  the  absence  of 
any  assumption  concerning  approximate  equidistribution  of  the  parameter  values  critical  for  F 
we  cannot  be  certain  of  the  number  of  bisection  steps  required. 

Megiddo  circumvents  this  difficulty  as  follows.  Suppose  that  a  can  be  executed  as  a 
parallel  algorithm,  which  given  an  input  array  A  (ji)  of  size  n  uses  0(p(n))  processors  in 
parallel,  and  that  together  these  perform  a  total  of  0{w(n))  elementary  computational  steps. 
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We  suppose  that  the  parallel  algorithm  alternates  between  phases  of  completely  asynchronous 
parallel  operation  and  brief  synchronization  cycles,  and  that  the  total  number  of  synchroniza- 
tion cycles  required  is  0(5(n)).  Suppose  that  the  elements  of  the  array  A  {\i)  are  used  only  in 
the  following  restricted  ways. 

(i)  They  are  combined  algebraically,  but  the  algebraic  combinations  B{\l)  that  are 
formed  are  of  uniformly  bounded  degree  in  terms  of  the  original  inputs  A  (p.). 

(ii)  The  quantities  5(jji)  which  result  are  compared  and  the  results  of  these  comparisons 
are  used  in  conditional  branches  by  the  processors  performing  the  comparisons.  However, 
no  processor  performs  more  than  0(1)  successive  comparisons  between  any  two  synchroniza- 
tion cycles  of  the  (parallel)  algorithm  a. 

In  this  case  we  can  proceed  as  follows.  Simulate  the  action  of  the  0(j>{n))  separate 
processors,  one  after  the  other.  Once  simulation  of  the  activity  of  a  particular  processor  has 
been  initiated,  it  continues  up  to  the  first  point  at  which  this  processor  either  compares  two 
quantities  fl(ji)  or  reaches  a  synchronization  point.  Algebraic  combination  of  inputs  A{\i.) 
and  intermediate  quantities  B{\x.)  to  produce  new  intermediate  quantities  are  performed,  but 
handled  symbolically.  As  soon  as  a  comparison  between  two  intermediate  quantities  fli(ji.), 
J}2(m-)  is  about  to  be  performed,  simulation  moves  on  to  the  next  processor  of  the  p  (n)  in 
turn,  until  the  comparisons  needed  to  be  performed  by  all  processors  have  been  collected. 
Since  both  the  quantities  Bi(^.)  and  BiiV')  involved  in  any  such  comparison  are  algebraic  in  \l 
and  of  degrees  not  exceeding  a  bound  known  in  advance,  the  points  p.],  ...,(!„  at  which 
any  single  comparison  like  B  i{\).)<B i{\i.)  can  reverse  are  bounded  in  number.  We  collect  all 
these  points,  as  they  are  generated  by  comparisons  attempted  by  successive  processors,  accu- 
mulating them  in  a  comprehensive  critical  point  list  Z,.  If  a  processor  performs  no  comparis- 
ons before  reaching  a  synchronization  point  of  the  parallel  algorithm,  simulation  of  its  activi- 
ties continues  until  the  synchronization  point  is  reached. 

Our  next  step  is  to  'resolve'  all  the  comparisons  which  this  simulation  procedure  has 
defined.    For  this,  we  sort  the  accumulated  list  L  of  critical  values,  find  its  median  JT,  and 
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evaluate  the  value  F  (ir)  =  a(A(jr))  in  normal  rather  than  simulated  fashion.  By  assumption, 
this  computation  tells  us  whether  the  ji*  we  seek  lies  above  or  below  the  value  jT.  If  the 
length  of  L  is  /,  a  sequence  of  log  /  such  evaluations  will  find  the  two  successive  jjl;,h/.i  in 
the  list  L  between  which  p,*  lies.  Note  that  since  the  quantities  B^V-)  being  compared  are  of 
algebraic  degrees  satisfying  a  uniform  a  priori  bound,  the  length  of  the  list  L  is  bounded  by  a 
constant  times  the  number  p  (n)  of  processors  employed.  The  manner  in  which  we  have 
defined  L  implies  that  every  one  of  the  attempted  comparisons  fl  i(fjL)<B2(M')  ^^^  ^  boolean 
value  -y  that  is  the  same  for  all  \i.  in  this  range;  the  value  7  can  plainly  be  determined  by  a 
purely  symbolic  calculation.  By  assigning  each  such  comparison  the  boolean  value  calculated 
for  it  in  this  way  we  resolve  all  required  comparisons  in  a  manner  appropriate  for  every  ji  in 
the  range  in  which  we  are  interested,  thus  making  it  possible  to  continue  our  simulation  until 
each  processor  either  attempts  one  more  comparison  or  advances  to  the  next  synchronization 
point. 

The  cost  of  this  simulation  is  the  sum  of  the  cost  O  (s  (n)w  *  (n)  log  p  (n))  of  all  required 
evaluations  a(A{\r)),  plus  the  cost  0{w{n))  of  the  'symbolic'  simulation  that  guides  the 
choice  of  a(A(fr))  to  evaluate.  The  expression  ^'(n)  occuring  in  the  cost 
0(s(n)w*(n)  log  pin))  denotes  the  cost  of  any  single  serial  or  parallel  evaluation  of  the  func- 
tion /"(p.),  since  each  such  evaluation  may  as  well  be  performed  serially.  The  (possibly  more 
expensive)  parallel  version  of  the  algorithm  is  required  only  to  guide  the  'symbolic'  simula- 
tion, in  particular  to  accumulate  the  critical  point  lists  L.  (Note  also  that  in  reckoning  the 
cost  H'(n)  of  the  parallel  algorithm  a,  we  can  omit  all  bookeeping  costs  incurred  merely  to 
assign  processors  to  tasks,  since  our  simulation  (which  is  serial)  will  not  incur  any  of  these 
costs.  That  is,  we  can  make  our  estimates  using  the  favorable  cost  model  of  parallel  compu- 
tation introduced  by  Valiant  [Va].) 

These  general  reflections  apply  in  the  following  way  to  the  case  that  concerns  us.  Con- 
sider a  single  pruning  phase  of  FINDMAX  (recall  that  it  has  O(log  n)  such  phases).  First 
group  the  planes  /,  into  disjoint  pairs,  and  for  each  such  pair  consider  the  dividing  line 
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L,j  =  (p,(.y^)-Pjiy-))ix,y,l)  =  0,  whose  coefficients  are  plainly  linear  functions  of  n.  The 
FINDMAX(n.)  algorithm  that  we  have  described  calculates  three  medians,  namely  the 
median  of  the  angles  formed  between  the  lines  and  the  positive  x-axis,  and  also  the  medians 
y„  and  *„.  Calculation  of  y„  is  typical,  since  the  other  two  required  median  calculations  are 
almost  identical  (from  the  point  of  view  of  our  generic  simulation).  It  is  easily  checked  that 
the  y  coordinates  of  the  points  of  intersection  of  the  pairs  of  lines  Ly  are  rational  quadratic 
functions  of  jt;  denote  these  functions  by  y\i\i),  ■  ■  .  ,y*(M-)-  To  apply  Megiddo's  technique, 
we  require  a  parallel  algorithm  for  finding  the  median  of  n  elements;  for  this  we  can  use  the 
recently  developed  algorithm  by  Ajtai,  Komlos,  Steiger  and  Szemeredi  [AKSS]  which  runs  in 
parallel  time  0(log  log  n)  and  uses  0{n)  processors,  or,  alternatively,  we  run  the  algorithm 
of  Cole  and  Yap  [CY],  which  runs  in  parallel  time  0((log  log  n)^)  and  uses  0(n)  processors. 
(The  constants  appearing  in  the  time  complexity  bound  of  the  [CY]  algorithm  arc  substan- 
tially smaller  than  those  of  the  [AKSS]  algorithm.)  These  algorithms  meet  the  conditions  for 
application  of  Megiddo's  technique;  e.g.  the  [AKSS]  algorithm  involves  O(log  log  n)  succes- 
sive synchronization  cycles,  between  which  each  of  the  0{n)  processors  used  performs  just  a 
single  comparison. 

If,  for  example,  we  apply  Megiddo's  transformation  to  the  [AKSS]  algorithm,  the 
resulting  procedure  finds  the  median  of  the  points  y,(M.)  in  the  correct  range  of  p,  (that  sur- 
rounding ji"^)  at  a  total  computational  cost  of  0{n  log  n  log  log  n).  The  other  median  calcu- 
lations required  proceed  in  essentisilly  the  same  manner. 

In  addition,  the  current  pruning  phase  involves  two  calls  to  LINETEST,  which  also  have 
to  be  simulated  generically.  Each  execution  of  LINETEST  involves  O(log  n)  pruning  phases, 
each  requiring  one  median  calculation,  followed  by  one  evaluation  of  F.  The  median  calcula- 
tion can  proceed  in  the  same  manner  outlined  above,  whereas  to  calculate  F  (i.e.  the 
minimum  of  the  functions  /j,(ji)(j:,y,  1)),  we  use  a  similar  approach  based  on  the  parallel 
minimum-finding  algorithm  of  Valiant  [Va],  which  also  uses  0(n)  processors  and  runs  in 
0(log  log  n)  parallel  time.  Thus  a  single  (generic)  execution  of  LINETEST  can  be  performed 
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in  time  0{n  log^n  log  log  n).  (This  bound  is  easily  seen  to  apply  also  to  all  the  remaining 
steps  of  LINETEST,  such  as  determining  the  side  of  the  line  containing  the  global  maximum 
of  F,  etc.) 

Since  the  execution  of  FINDMAX  involves  0(log  n)  pruning  stages,  we  conclude  that 
the  generic  execution  of  FIND-ji  can  be  performed  in  overall  time  0{n  log^n  log  log  n). 

When  FIND-p.  terminates,  we  will  have  found  a  range  for  jji*  with  the  property  that  the 
sequence  of  steps  that  the  maximum  finding  procedure  FINDMAX  would  go  through  if 
applied  to  any  p.  in  this  range  would  be  independent  of  ^l..  We  can  then  take  p,'^  to  be  the 
smallest  endpoint  of  this  range. 

Calculation  of  ji"  proceeds  in  fully  symmetric  manner;  the  maximum  of  ji*  and  ji"  is 
the  value  of  the  coefficient  of  friction  that  we  seek.  This  completes  the  description  of  our 
algorithm.    Its  computational  cost  is  plainly  0(n  log^n  log  log  n). 

Remark:  The  maximum-finding  procedure  can  readily  be  modified  to  produce  (one  or  two 
or)  three  planes  along  whose  intersection  the  maximum  of  M  is  attained.  Hence  it  can  also  be 
used  to  represent  the  origin  as  a  convex  combination  of  (two  or)  three  of  the  points  p,(M-*)- 

We  summarize  the  results  of  this  section  in  the  following 

Theorem  2,1:  Given  n  fingers  grasping  a  rigid  free  2-D  object  B  at  specified  points  with 
specified  normals,  the  smallest  coefficient  of  friction  p,*  under  which  a  stable  grasping  of  B 
by  the  fingers  is  possible,  as  well  as  the  force  targets  to  be  exerted  by  the  fingers  to  realize 
such  a  grasp,  can  be  computed  in  time  0(n  log^n  log  log  n). 

Remark:  We  can  slightly  improve  the  running  time  of  the  algorithm  using  the  recent  tech- 
nique of  Cole  [Co]  for  speeding-up  Megiddo's  method.  However,  Cole's  technique  cannot  be 
applied  in  Valiant's  model  of  computation,  since  it  requires  the  underlying  parallel  algorithm 
to  be  implementable  on  a  network  (i.e.  with  prespecified  I/O  coimections,  of  low  in-degree, 
between  processors  in  adjacent  parallel  steps  of  the  algorithm).  Using  the  [AKS]  sorting  net- 
work rather  than  the  median  or  minimum  finding  procedures  that  we  have  used  above,  one 


-  17- 

can  show  that  Cole's  technique  improves  the  running  time  of  the  algorithm  down  to 
0(n  log^n). 

3.    Calculating  finger  force  targets  for  a  given  external  force  and  torque 

The  procedure  FIND-jt  described  above  finds  the  'optimal'  frictional  grip  that  can  be 
exerted  on  a  body  not  subject  to  outside  forces.  However,  in  many  applications,  gripping 
robot  fingers  will  need  to  manipulate  a  body  subject  to  varying  external  forces  and  torques 
(e.g.  in  the  presence  of  gravity,  or  when  the  body  makes  contact  with  another  object.) 
Assume  that  some  sufficiently  accurate  measurement  of  these  external  forces  and  torques 
(Fe,Te)  is  available  at  each  moment.  A  reasonable  control  heuristic  is  then  to  find  the  smallest 
coefficient  of  friction  p.*  =  p.*  for  which  the  3-vector  e  =  (f  ,,t,)  can  be  expressed  as  a  posi- 
tive combination  of  the  points  p  i(m-').  •  ■  •  .?«(>»•*).  ie.  the  smallest  jl*  for  which  e  belongs 
to  the  cone  C(m.').  Since  at  this  ji*  the  vector  e  will  lie  on  the  boundary  of  C(ji,*),  hence  on 
a  face  of  that  cone  generated  by  the  origin  and  two  points  p,(m.*),  Py(^.*),  it  follows  that  in 
general  e  will  be  expressible  as  a  positive  combination  of  just  two  of  the  vectors  p/  at  the 
optimal  jA* . 

Remark:  Unlike  the  case  studied  in  the  previous  section,  here  p.,  is  not  always  well  defined, 
namely  the  infimum  jji*  of  all  values  of  p,  for  which  e  €  C(p.)  need  not  satisfy  this  condition. 
Indeed,  as  n  increases  and  C  (ti)  expands,  there  will  be  a  critical  level  ti,^  at  which  C(ji.)  is  a 
half-space  bounded  by  a  plane  passing  through  the  origin.  If  e  is  any  vector  in  the  open  com- 
plementary half-space,  then  e  €  C(m.)  for  any  ^>.>^x.*  but  e  is  not  in  C(fi,*).  However,  for 
each  e  €  C{y.'*^),  the  quantity  ji,  is  well  defined.  This  observation  is  reflected  in  the  subse- 
quent analysis. 

For  each  jj,  <  ji*  the  vector  e  lies  outside  the  cone  C(p,).  It  is  easy  to  check  that  this  is 
equivalent  to  the  existence  of  a  plane  ir  separating  e  from  C(ji.)  and  supporting  C(ji)  at  the 
origin.  Let  x  be  the  vector  normal  to  it  and  pointing  into  the  side  of  ir  containing  C{^).). 
Then  we  must  have 


18  - 


ex  s  0 
p,(n)x  S:  0  ,         i  =  I,  .  .  .  ,n  . 
Two  approaches  to  the  problem  suggest  themselves.  One  is  to  calculate 

n-  1 

max   min  p,(n)x  , 

«       1  =  1 

where  we  have  put  p„^i(ji)  =  -e,  and  then  find  the  smallest  value  of  ft  (or  rather  the 
infimum  of  all  values  of  ji)  for  which  this  maximum  is  nonpositive.  This  can  be  done  in 
exactly  the  same  way  as  in  the  preceding  section,  at  computational  cost  0(n  log^/i  log  log  n) 
(or  0{n  log'n)).  However,  this  solution  requires  application  of  FIND-ji  each  time  the  exter- 
nal force-and-torque  vector  e  changes,  i.e.  once  each  control  cycle,  and  is  therefore  prohibi- 
tive (even  from  a  theoretical  point  of  view).  Alternative  approaches  are  therefore  desirable. 
These  can  involve  (possibly  costly)  preprocessing  before  object  manipulation  begins,  and  then 
use  the  structures  generated  by  preprocessing  to  find  the  optimal  ji.  and  the  balancing  finger- 
force  targets  quickly  during  each  control  cycle. 

This  can  be  done  as  follows.  Consider  the  family  of  convex  cones  C(jjl)  spanned  by  the 
vectors  p,(ji).  These  increase  steadily  with  jx,  from  the  minimum  cone  C(0).  As  long  as 
C(fji)  is  not  all  of  3-space  it  must  lie  entirely  on  one  side  of  some  plane.  As  in  the  remark 
made  above,  let  n*  be  the  smallest  value  above  which  C(ji)  is  the  whole  of  3-space,  let  Pq 
be  a  plane  through  the  origin  such  that  C(jjl*)  lies  entirely  on  one  side  of  Pq.  and  let  n  be  the 
vector  normal  to  Pq  and  pointing  toward  the  side  H  of  Pq  containing  C(jji'^).  Then  plainly 
nxaO  for  all  x  €  C(m.),  ^.s^.*.  Take  any  fixed  plane  P  parallel  to  Pq  and  intersecting 
C(ii).  Let  e  be  any  vector  (of  external  force/torque)  pointing  into  the  halfspace  H,  and  let 
e  €  P  be  an  appropriate  constant  multiple  of  e.  Plainly  e  €  C{\l)  if  and  only  if  e  €  C(m.). 
Thus  our  problem  is  to  find  the  smallest  p.  for  which  e  i  C(n).  For  each  /  let  9,(m.)  be  the 
constant  multiple  of  p,(n)  which  lies  in  P.  Then  e  €  C(n)  iff  e  belongs  to  the  convex  hull 
K^ii)  of  the  points  9,(jji.);  thus  our  problem  is  to  determine,  for  each  point  e  of  the  plane  P, 
the  smallest  ^  for  which  e  belongs  to  K(it,).  (Note  that  this  set-up  makes  sense  only  when  e 
lies  in  H.    However,  if  e  lies  in  the  other  side  of  Pq  we  can  take  p.*  to  be  any  number  greater 
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than  ii.^  (in  the  spirit  of  the  remark  made  above)  and  find  four  fixed  vectors  Pi(m.*), 
P2(m-*).  PsCM'*).  P4(M'*).  whose  cone  spans  the  entire  3-space,  thus  allowing  us  to  express 
any  such  e  as  a  positive  combination  of  two  or  three  of  these  vectors.) 

The  convex  polygons  Ar(ji.)  increase  steadily.  Moreover,  since  the  coordinates  of  the 
points  p,(jji)  depend  linearly  on  p.  the  coordinates  of  their  'projections'  q,(yi)  are  rational 
linear  functions  of  ji.  For  each  ^l  some  subset  of  the  points  9,(ji.)  will  comprise  the  vertices 
of  K(it.).  Call  each  value  of  ji  at  which  this  set  changes  a  critical  level  of  p..  Clearly  at  each 
such  critical  level  either  two  of  the  vectors  p,(m-)  must  be  collinear  or  three  must  be  coplanar; 
since  in  our  application  the  pXi'-)  depend  linearly  on  p.,  these  conditions  can  be  expressed  by 
algebraic  equations  of  degree  at  most  three,  so  the  number  of  critical  levels  is  at  most  O(n^). 
In  fact,  as  will  be  shown  below,  the  actual  number  of  these  levels  is  only  little  more  than 
quadratic  in  n. 

Between  any  two  successive  critical  levels  ji(<p.^  a  fixed  set  of  the  points  ?,(j».)  will 
comprise  the  vertices  of  Ar(ji.),  and  these  will  occur  in  some  fixed  circular  order.  The  straight 
segments  traversed  by  the  vertices  9,(ji.)  divide  the  polygonal  annuli  K(\i.>,)  -  K{y.i)  into  at 
most  0{n)  regions  /?,_,  each  bounded  by  two  successive  such  segments.  Each  e  in  R,j  can 
plainly  be  written  as  a  convex  combination  e=a^,(^^.)+ (l-a)^^(jjL)  for  some  ji  between  jjl^ 
and  ^/,  and  (one  of)  the  p.  appearing  in  a  representation  of  this  kind  is  the  smallest  p.  for 
which  e  can  be  written  as  a  positive  combination  of  the  vectors  Pi(n.)-  The  required  \i  can  be 
easily  computed  by  noting  that  the  vectors  e,p,(jji),pj(jji)  must  be  coplanar  for  that  p.,  which 
implies  that  the  relevant  p,  can  be  found  by  solving  the  quadratic  equation 
det  (e,p,(fJt).P;(p^))  =  0. 

Accordingly,  we  can  construct  an  algorithm  of  the  following  form.  Build  an  ordered  list 
L  of  all  the  critical  levels  n,i,ji,2,  .  .  .  ,  of  p..  Then  partition  P  into  disjoint  polygonal  regions 
by  drawing  the  'critical'  convex  sets  ^(jjl,),^((12),  .  .  .  ,  and  by  further  partitioning  each 
annulus  between  two  successive  A^(ji,,)'s  into  the  polygonal  regions  R,j  defined  above  (note 
that  each  /fy  is  cither  a  triangle  or  a  quadrilateral).    Build  up  the  data  structure  required  to 
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support  an  O(log  n)  point  location  algorithm  for  determining  the  region  R,j  in  which  any 
given  e  €  P  lies  (using  e.g.  the  techniques  of  [Ki],  [EGS]  or  [ST]).  Given  any  e,  apply  this 
point  location  algorithm  to  determine  i,j  and  then  solve  the  quadratic  equation 
det  (e,Pi(jj,),p;(M.))  =  0  to  determine  the  desired  jji  =  jie. 

The  space  requirements  and  preprocessing  cost  of  this  procedure  depend  on  the  number 
K(n)  of  critical  levels  of  jt  at  which  the  convex  hull  of  the  n  2-D  vectors  9i(m.),  .  .  .  ,q„iii.) 
changes  combinatorially.  This  quantity  was  estimated  by  Atallah  [At]  who  showed  that  the 
number  of  such  levels  is  0(nX.j(n)),  where  Xj(n)  is  the  maximum  length  of  an  {n,s)  Daven- 
port Schinzel  sequence,  which  is  also  equal  to  the  maximum  number  of  graph  portions  which 
can  appear  in  the  lower  envelope  of  n  continuous  functions,  each  pair  of  which  is  allowed  to 
intersect  in  at  most  5  points.  (The  integer  s  that  applies  in  our  case  will  be  characterized 
below.)  It  is  known  that 

(i)      \,(n)  =  0(n)  for  J  =  1,2, 

(ii)     X3(n)  =  6(na(n)),  where  a  is  the  very  slowly  growing  inverse  of  the  Ackermann  func- 
tion [HS], 

(iii)  X4(n)  =  e(n-2'»(">)  [ASS], 

(iv)  \2,(n)  =  0(n-2^(»"'>'"^>),  for  s>2  [ASS], 

(v)  X2,(n)  =  n(n-2"<»(">'"^>),  for  s>2  [ASS], 

(vi)  Xz.  +  iC")  =  0(n-a(n)°(»<'"'"'>),  for  s^2  [ASS]. 

In  the  following  pages  we  will  use  a  very  similar  method  to  derive  an  upper  bound  for  k(r), 
together  with  an  algorithm  for  the  calculation  of  the  critical  levels  )i,  and  the  associated  parti- 
tioning of  P.  The  technical  differences  between  our  discussion  and  that  of  [At]  are  required 
to  handle  various  difficulties  arising  from  the  manner  in  which  the  points  9^((i)  depend  on  (t 
in  our  case. 

The  estimates  we  require  can  be  derived  as  follows.  Since  the  sets  K(ii)  expand  with 
ii,  we  can  assume  without  loss  of  generality  that  all  these  sets  contain  the  origin.    For  each 
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Q  g  (_  JL  .21)  and  each  ji,  let  F,(6,jji.)  be  the  distance  from  the  origin  to  the  line  passing 

through  the  i-th  point  9,(fi)  and  perpendicular  to  the  vector  (cos  6,  sin  6).  More  precisely, 
put 

^,(e,f«.)  =  iq,M)x  cos  e  +  iq,M)y  sin  6  . 
Let 

A/(e,ji)  =  maxF,(e,ji) 
be  the  upper  envelope  of  the  functions  f,  thereby  defined.  Then  clearly  for  each  e.ji.  the 
value  A/(e,M.)  is  the  distance  from  the  origin  to  the  line  supporting  ATCm-)  and  normal  to 
(cos  6,  sin  6).  Thus  calculation  of  the  convex  hulls  K  (ii)  is  equivalent  to  the  calculation  of  M; 
the  critical  values  of  ji  are  precisely  those  p.  for  which  the  combinatorial  structure  of  the 
upper  envelope  function  Af  (♦,fj.)  changes. 

It  is  easily  seen  that  a  point  (6,^.)  at  which  three  of  the  vertices  of  JC^(ji)  become  col- 
linear  is  also  a  point  at  which  the  maximum  M(0,(jl)  is  attained  simultaneously  by  three  (or 
more)  functions  F,,  Fj,  Fk,  and  conversely.  Let  (9,>i.)  be  such  a  point.  Then  for  one  of  the 
three  (or  more)  indices  i,  j,  k,  which  without  loss  of  generality  we  can  take  to  be  /,  we  must 
have  F,(ji,,e')  >  Fj(n,,Q')  and  /".(jji.e')  >  /'^(ji.e')  for  each  6'  larger  than  6  and  sufficiently 
near  it. 

For  each  pair  ij  consider  the  curve  7,;  in  the  6, ft  plane  at  which  F,  =  Fj.  For  each 
fixed  (ji  this  condition  defines  a  unique  value  6  =  ^^(m.)  determined  by  the  intersection  of 
two  straight  lines.  For  angles  6'  somewhat  to  the  left  of  this  intersection  point  we  have  con- 
sistently either  F,(e',ji.)  >  F;(e',ji)  or  the  reverse,  and  vice  versa  for  9'  to  the  right  of  the 
intersection  point. 

G^"(u) 
Since  the  function  -y^Xji)  can  also  be  written  as  9  =  arctg  —^ where  2^'\  Q'-^^  are 

Gy  (m-) 
quadratic  functions  of  ^,,  it  has  at  most  two  points  of  (jump)  discontinuity  and  is  defined  and 
continuous  at  all  other  points.   In  what  follows  we  will  consider  each  connected  portion  of  the 
graph  of  y,j  as  a  separate  (partially  defined)  function  of  ^l.,  which  still  leaves  us  with  at  most 
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0(n)  such  functions  for  each  fixed  i. 

Suppose  6  =  "iijiti)  is  defined  for  a  particular  pair  i,j.  Then,  as  remarked  above, 
throughout  the  entire  domain  of  definition  /  =  I,j  of  the  (partial)  function  y,j,  we  must  have 
either 

(i)  f,(e,fi.)  >  Fj(B,yi)  (resp.  F,(e,^.)  <  F;(e,M.))  for  each  jt  €  /  and  each  0  >  7,j(ft)  (resp. 

0  <  7y(^l)).  or  else 

(ii)  F,(e,ji.)  <  Fj(6,ii.)  (resp.  F,(e,M.)  >  f;(e,n))  for  each  p,  €  /  and  each  6  >  y,j(y.)  (resp. 

e  <  y.jM)- 

Consider  a  particular  ji,  and  suppose  without  loss  of  generality  that  we  deal  with  the 
first  case.  Then  clearly  for  each  ji.  €  /  and  each  6  <  "y,;(fi.),  the  function  F,  lies  below  the 
maximum  value  Af  (6,^,)  =  max  ^^(e.ji)  at  the  point  (e.jji).  Hence  if  we  put 

«i(,(M.)  =  max  7,;(fi) 
where  the  maximum  is  taken  over  all  j  for  which  p.  €  7,^  and  F,{B,^^,)  >  Fy(e,|i)  for  each 
6  >  ■Yij(m-).  it  follows  that  F,  can  only  attain  the  maximum  value  M(8,ji)  at  a  given  (6,fi.)  if 
e  2:  i|;,(ji.). 

We  noted  above  that  at  any  point  (e.ji.)  at  which  three  or  more  of  the  vertices  of  K{\i) 
become  collinear  there  exists  an  i  such  that  F,(6,n)  =  Af  (O.ji)  while  F,(9',ji,)  >  Fj{Q' ,\jl) 
for  ally  7^1  and  e'>6  which  is  sufficiently  near  6.  Moreover  there  exist  at  least  two  indices  y 
and  *  distinct  from  i  such  thatF,(e,M.)  =  Fj(9,ii)  =  Fi(9,n). 

From  this  it  is  easy  to  see  that 

0  =  luM  =  7<*(M')  =  <I'.(M') 
for  at  least  two  distinct  j  and  k  simultaneously.    Hence  every  p,  for  which  three  vertices  of 

the  convex  hull  become  collinear  must  be  a  'breakpoint'  at  which  at  least  one  of  the  maximiz- 
ing 'envelope'  functions  <|»,(ji.)  =  max  7,;(m.)  is  simultaneously  the  value  of  at  least  two  func- 
tions yij,  -y,^.  Thus  calculation  of  the  collection  of  all  breakpoints  in  the  graphs  of  all  the  i|i, 
gives  us  a  superset  of  the  set  of  all  critical  levels  of  ji  at  which  three  vertices  of  K(^J.)  become 
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collinear. 

The  combinatorial  results  of  [At],  [HS],  [ASS]  allow  us  to  estimate  the  number  of  these 
'breakpoints'  by  noting  that  any  two  curves  -y,;,  -y,*  can  intersect  in  at  most  three  points. 
Indeed,  since  each  point  of  intersection  of  these  curves  is  a  point  at  which  the  three  functions 
F,,  Fj,  Fk  must  coincide,  and  hence  a  p.  for  which  the  three  points  p,{ii.),  PjM,  Pt(M-)  must 
be  collinear.  Hence  the  parameter  values  ji.  for  which  such  an  intersection  can  occur  satisfy  a 
cubic  equation.  It  follows  that  the  two  curves  y,j,  y,k  can  intersect  in  at  most  three  points.  If 
the  functions  -y/;  were  total  and  continuous,  it  would  follow  at  once  that  the  number  of  break- 
points in  i|i,  would  be  at  most  XsCn),  as  defined  above.  However,  since  each  of  the  partial 
functions  lij  with  which  we  need  to  work  may  only  be  defined  on  a  subinterval  of  the  real 
axis,  this  argument  only  allows  us  to  bound  the  number  of  breakpoints  by  the  somewhat 
larger  (but  stiU  nearly  linear)  quantity  XjCn)  =  0(na(n)^<'"^''»)  (cf.  [At],  [ASS],  [SS]  for 
details).  Overall  the  total  number  of  critical  p.  levels  is  at  most  0(n\s(n)),  i.e.  little  more 
than  quadratic  in  n. 

As  described  in  [At],  [HS],  these  breakpoints  in  each  of  the  maximizing  'envelope'  func- 
tions «|«,  can  be  calculated  using  a  divide  and  conquer  approach.  Roughly  speeiking,  for  each  i 
we  first  collect  all  functions  -y/y  which  satisfy  condition  (i)  above.  Then  we  divide  the  collec- 
tion of  these  functions  into  two  roughly  equal  subsets,  recursively  calculate  the  maximum  of 
the  functions  in  each  subset  separately,  and  then  merge  these  two  maxima  in  the  straightfor- 
ward manner  described  in  [At],  [HS].  The  complexity  of  this  procedure,  summed  over  all 
envelopes  \|i,,  is  bounded  by  0{n\i{n)  log  n),  and  as  already  noted  the  procedure  yields  a 
superset  of  all  critical  values  of  p,. 

For  full  calculation  of  the  smallest  ji  for  which  e  €  C  (jt)  we  also  need  to  calculate  the 
polygonal  partitioning  of  the  plane  P  described  above,  or  something  equivalent  to  it  for  our 
purposes.  It  is  more  efficient  to  work  with  the  slightly  modified  partitioning  of  P  defined  as 
follows.  Suppose  that  for  some  ji  the  convex  hull  K{^i.)  has  ?/(p.),  9;(m.)  as  two  successive 
vertices.  Let  /  =  (p,;,n;,)  be  the  maximal  interval  containing  p,  for  which  this  property  holds. 
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(Clearly  /  is  delimited  by  two  critical  levels  of  jx,  but  these  two  levels  need  not  be  adjacent  in 
the  list  of  critical  levels.)  Partition  P  into  disjoint  quadrilaterals  Q,  each  of  which  is  defined 
by  two  points  9,(fi.),  qj{\i.)  and  by  an  interval  /  =  (ji/.jn)  of  the  kind  just  described.  Thus  Q 
is  bounded  by  the  two  segments  connecting  (?,(ji/)  to  9;(n;)  and  qXv^h)  to  9y(M-A)  respectively 
and  by  the  two  straight  orbits  of  qX\i.),  qj{,\i.)  traversed  as  \i.  varies  through  /.  Since  (under 
the  technical  assumption  that  no  more  than  three  vertices  of  K{\l)  become  coUinear  at  any 
critical  level)  each  critical  level  of  p.  will  delimit  at  most  three  such  intervals  (associated  with 
the  three  pairs  of  points  9,  which  become  collinear  at  this  level),  it  follows  that  the  number  of 
these  quadrilaterals  is  at  most  0(nX5(n)). 

These  regions  can  be  calculated  as  follows.  First  build  a  sorted  list  L  of  (the  previously 
described  superset  of)  all  critical  levels  of  \i..  Process  these  critical  levels  in  increasing  order, 
maintaining  a  list  of  'active'  regions  whose  two  'lower'  vertices  q,(,\x.i)  and  qj{\i.i)  have  already 
been  encountered,  but  whose  two  remaining  vertices  q,{V'h)  and  qj{\ih)  have  not  yet  been 
encountered.  At  each  critical  level  pn,  we  'close'  the  small  number  of  active  regions  affected 
by  the  vertex  coUinearity  that  occurs  at  n*,  and  make  active  a  few  new  quadrilateral  regions 
of  the  same  kind.  This  procedure  generates  the  desired  partitioning  of  P,  and  plainly  can  be 
implemented  to  run  in  time  0{n\(,{n)  log  n). 

Once  the  quadrilateral  partitioning  generated  by  this  procedure  is  available,  we  can  use 
any  one  of  the  fast  point  location  preprocessing  algorithms  mentioned  above  to  set  up  the 
required  data  structures  for  logarithmic-time  location  queries. 

In  summary,  we  have 

Theorem  3.1:  The  0(n)  finger  force/torque  vectors  p,(|ji)  can  be  preprocessed  in  time 
0{n\i{n)  log  n)  to  produce  a  data  structure  of  size  0(nX5(n)),  from  which,  in  C>(log  n) 
time,  we  can  determine  the  smallest  coefficient  of  friction  \x.  for  which  any  given  external 
force/torque  vector  e  is  expressible  as  a  positive  combination  of  (two  of)  the  vectors  p,{,ii.), 
and  also  produce  such  a  positive  combination  when  it  exists. 
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4.   Pragmatic  techniques  applicable  when  only  a  small  number  of  Hngers  are  involved. 

We  close  the  paper  with  a  few  observations  concerning  practical  solutions  to  the  prob- 
lems studied  above.  Even  though  the  solutions  that  we  have  developed  in  the  preceding  sec- 
tions are  mainly  of  theoretical  interest,  the  geometric  analysis  that  accompanies  them  can  be 
used  to  obtain  variants  of  these  algorithms  which  are  useful  in  practice. 

Consider  first  the  problem  of  balancing  a  free  body,  as  studied  in  Section  2.  Our 
analysis  suggests  the  following  simple  way  to  find  ft*.  Calculate  all  critical  levels  of  ji  at 
which  three  vectors  p.Cm-)  become  coplanar,  and  perform  a  binary  search  through  them  to 
find  p,*,  using  the  linear  programming  based  strategy  of  Section  2  (or  perhaps  some  simpli- 
fied variant  of  it)  at  each  binary  search  step.  Since  there  are  O(n^)  critical  levels,  this  solu- 
tion will  run  in  0(n^log  n)  time,  which,  when  n  is  small,  can  be  quite  reasonable.  An  alter- 
native method  is  simply  to  perform  binary  search  over  p.,  halving  the  range  for  n,*  at  each 
step,  until  ji*  is  determined  to  sufficient  accuracy. 

Next  consider  the  problem  studied  in  Section  3.  The  solution  presented  there  is  not  too 
far  from  being  practical.  First,  it  provides  a  general  framework  for  solving  the  problem. 
Second,  the  analysis  given  there  is  useful  in  establishing  that  the  number  of  critical  levels  of 
(1.  is  relatively  small.  To  obtain  a  practical  algorithm  out  of  the  one  given  in  Section  3,  we  can 
simply  calculate  the  polygonal  subdivision  described  there,  and  store  it  in  any  convenient 
manner  which  facilitates  fast  point  location  in  practice,  using  less  fancy,  but  easier  to  manipu- 
late, data-structures.  This  approach  seems  to  us  to  be  ideal  for  situations  that  require 
repeated  computations,  within  a  control  loop,  of  force  targets  for  the  fingers  needed  to  bal- 
ance varying  external  force  and  torque. 
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